Sticky steps inhibit step motions near equilibrium 
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Using a Monte Carlo method on a lattice model of a vicinal surface with short-range step-step 
attraction, we show that, at low temperature and near equilibrium, there is an inhibition of the 
motion of macro-steps. This inhibition leads to a pinning of steps without defects, adsorbates, 
or impurities (self- pinning of steps). We show that this inhibition of the macro-step motion is 
caused by faceted steps, which are macro-steps that have a smooth side surface. The faceted steps 
result from discontinuities in the anisotropic surface tension (the surface free energy per area). The 
discontinuities are brought into the surface tension by the short-range step-step attraction. The 
short-range step-step attraction also originates 'step-droplets', which are locally merged steps, at 
higher temperatures. We derive an analytic equation of the surface stiffness tensor for the vicinal 
surface around the (001) surface. Using the surface stiffness tensor, we show that step-droplets 
roughen the vicinal surface. Contrary to what we expected, the step-droplets slow down the step 
velocity due to the diminishment of kinks in the merged steps (smoothing of the merged steps). 

PACS numbers: 81.10.Aj, 68.35.Ct, 05.70.Np, 68.35.Md 



I. INTRODUCTION 

A vicinal surface is a tilted surface deviated from a 
low Miller-index surface. The vicinal surface is thought 
to be described by the terrace-step-kink (TSK) modelpj- 
Q when the temperature is below the roughening transi- 
tion temperature of the low Miller-index surface. When 
the surface steps are regarded as linear excitations, the 
surface free energy per projected x-y area f(p)[&, 7] is ob- 
tained as the ground-state energy of the one- dimensional 
(ID) free fermion (FF) in the following form[l|-[16|: 




f(p) = f(0)+fP+Bp 3 + O(p% 



(1) 



where p represents the step density, 7 represents the step 
tension, and B represents the step-interaction coefficient. 
The characteristic feature of Eq. ([IJ is the lack of a 
quadratic term for p. The form of the free energy of Eq. 
([T]) is common to other many-body systems with nonover- 
lapping linear excitations embedded in two dimensions 
(2D) [8- 10]. The p-expanded form of the free energy is 
called the Gruber-Mullins-Pokrovsky-TalapovJGMPT) 
universal formQ or the ID FF universal form (l2i. [Til [l5j . 

If we introduce the surface gradient p = (p Xl p y ) = 
dz(x,y)/dx,dz(x,y)/dy), p is described by p as p = 
\p\/di, where d\ (=1) is the step height of an elemental 
step, and z(x,y) represents the surface height. We call 
f(p) the vicinal surface free energy. Using the surface 
gradient, the surface tension 7 sur f(n) links to the vicinal 
surface free energy as 



FIG. 1: (a) Perspective view of a surface on the RSOS model 
tilted towards the [100] direction, (b) Top view of a surface 
on the RSOS model tilted towards the [110] direction. Gray 
lines: surface steps. Filled squares: Collision points of the 
adjacent steps. 



7surf(n) = m/y/T+W, ( 2 ) 
where n is the normal unit vector of the surface at z(x, y). 
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In our previous papers [17H19j . we presented a sim- 
ple lattice model for a surface with sticky steps, the re- 
stricted solid-on-solid model with the point-contact-type 
step-step attraction (p-RSOS model) (Fig. [T]). Using 
methods of statistical mechanics, we determined the sur- 
face free energy with full anisotropy for the vicinal sur- 
face around the (001) surface. At low temperature, we 
obtained the anisotropic surface free energy, which in- 
cluded the discontinuities^^ due to the step-step at- 
traction. In addition, in one temperature region, we ob- 
tained non-GMPT values of the shape exponents fl9| on 
the profile of the equilibrium crystal shape (ECS), which 
is the shape of a crystal particle that has the least surface 
free energy 0-@],[2(J-[2i|. The non-GMPT values of the 
shape exponents obtained by the methods of statistical 
mechanics closely relate to the non-GMPT \p\ expanded 
form of the vicinal surface free energy, as follows 

fe S {f) = /(O)+7(^)b1+^c ff (^)bT + Bc ff (0)|p1 3 +O(bT). 

(3) 

Here we have introduced a polar coordinate system with 
respect to p, and we have p x = \p\ cos <j> and p y = \p\ sin <f>. 
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The angle <j> describes the angle of tilt from the y-axis, 
which is formed by the mean running direction of the 
steps. As seen from Eq. (j3|), a quadratic term with re- 
spect to |p| has appeared. In order to obtain Eq. ([3]), 
we took the inhomogeneity of the step configuration, the 
' step- droplets 1 (ID boson n-mers with finite lifetimes), 
into consideration]!!)]. The step-droplets are the locally 
merged steps in the vicinal surface. Using Eq. ([3]), we 
were able to reproduce the singularity in the surface free 
energy and the non-GMPT shape exponents, which were 
obtained with methods of statistical mechanics. 

Though there are extensive studies on the mechanisms 
of step bunchingfsuch as electromigration|2 7|- [29l| . the 
Schwoebel effectM |30| , the Shockwave effect [3 ll] Timpu- 
rity effects [Hj]- (34[, and strain effects 3 5j-[4]|). the vic- 
inal surface free energy is assumed to have the GMPT 
form (Eq. ([1}). A discontinuity in the anisotropic sur- 
face tension is expected to strongly affect the morphology 
of the crystal surface, accompanied by a change in the dy- 
namics of surface steps. The combination of the above 
mentioned mechanisms and a discontinuity in the surface 
tension will give us a true understanding of the nano- and 
mesoscale phenomena on the surface. In addition, since 
the self-assembling of surface steps is expected to provide 
new technology to construct nanoscale structures (such 
as nanowire) , the study of the effects of the discontinuity 
in the surface tension to the morphology and the surface 
dynamics is important from an industrial point of view. 

The aim of this paper, therefore, is to study the effect 
of the discontinuity in the surface tension to the step dy- 
namics near equilibrium. We also discuss the morpholog- 
ical aspect of the crystal surface from the viewpoint of the 
self- assembling of steps. In order to obtain clear results 
for the effect of the discontinuity in the surface tension, 
we exclude realistic interactions on the vicinal surface 
except for the point-contact-type step-step attraction. 

This paper is organized as follows. In §11, we present 
the restricted solid-on-solid (RSOS) lattice model with 
point-contact-type step-step attraction (the p-RSOS 
model) . Using a Monte Carlo method with the Metropo- 
lis algorithm, we demonstrate the inhibition and the 
slowing down of the step motion near equilibrium on the 
p-RSOS model. In §111, we show the discontinuous vici- 
nal surface free energy and the equilibrium configuration 
of the vicinal surface. We also calculate analytically the 
surface stiffness tensor near the (001) surface. In §IV, the 
inhibition of the motion of a macro-step, the intermittent 
motion of the vicinal surface, and the pinning phenom- 
ena in the vicinal surface are studied in connection with 
'step faceting' [42j. In §V, we study the roughness of the 
vicinal surface with step-droplets and the slowing down 
of the step motion near equilibrium. Summary and fur- 
ther discussion are given in §VI, and the conclusion is 
given in §VTI. 



II. SURFACE MOTIONS NEAR EQUILIBRIUM 

A. A lattice model for sticky steps 

In this section, we show Monte Carlo results for the 
surface motion near equilibrium for a vicinal surface 
with sticky steps. The model we adopted for the Monte 
Carlo calculation is the restricted solid-on-solid (RSOS) 
model with point-contact-type step-step attraction (p- 
RSOS model) 

To describe microscopic surface undulations, let us 
consider the surface height h(i,j) at a site (i,J) on a 
square lattice (Fig. [[]). In the RSOS model[H], the 
height differences between nearest-neighbor (nn) sites are 
restricted to values of {1,0, —I}. We consider a point- 
contact-type microscopic step-step interaction and refer 
to this model as the p-RSOS mo del [HI] -[II]. The Hamil- 
tonian for the p-RSOS model can then be written as 

^p-rsos = '^24\K i + - h (i>j)\ 

i,j 

+\h(i,j + l)-h(ij)\} 

+ J2e iat [S(\h(i + l,j+l)-h(i,j)\,2) 

+S(\h(i + l,j-l)-h(i,j)\,2)}, (4) 

where e represents the microscopic ledge energy, e ln t is 
the microscopic step-step interaction energy, and S(a,b) 
is the Kronecker delta. In the case of e; n t < 0, the in- 
teraction between the steps becomes attractive [8(| . The 
summation with respect to is performed over all 

sites on the square lattice. The RSOS restriction is re- 
quired implicitly. 

Physically, we consider that the point-contact-type 
step-step attraction arises from the local formation of a 
bonding state at the collision point of the adjacent steps. 
When adjacent steps collide at a point, the orbital of the 
dangling bond of each step will overlap, and the spin pair- 
ing between electrons in the dangling bonds will form a 
bonding state (45j. The energy gain obtained by forming 
the bonding state amounts to — e m t- 

In the p-RSOS model, there are two characteristic tem- 
peratures Tf^i and TV 2 for the vicinal surface tilted to- 
wards the [110] direction (4> = 7r/4) 19]. At tempera- 
tures T < Tf,i, the surface tension of the (111) surface 
becomes discontinuous. At temperatures T < T/.2, the 
surface tension of the (001) surface becomes discontinu- 
ous. That is, on the profile of the ECS, the surface slope 
of the vicinal surface jumps at the (111) facet edge (a 
first-order shape transition) for T < Tf t i, and the sur- 
face slope jumps at the (001) facet edge for T < T/.2. 
Then, the (001) facet directly contacts the (111) facet 
for T < T fi2 . 

The key concept that is introduced to explain the 
breakdown of the ID FF picture is the inhomogeneity 
resulting from the discontinuity in the surface tension. 
In other words, the formation of a 'step-droplet' (a ID 
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boson n-mer), which is formed from locally merged steps. 
Using the size of the step-droplet, A cS ((j)) and B e s((j>) in 
Eq. (|3]) are expressed explicitly in the limit of \p\ —> as 
follows fill: 



(5a) 



1 



2di 



(5b) 



where rfi (= 1) represents the height of an elementary 
step, and -Bi represents the step-interaction coefficient 
between the elementary steps. The variables Uq (4>), 
~f[ m \<fi), and j[ m \<p) are defined by 



a m ( 7 „(^)/n) 



dn rn 
d m {ln{4>)/n) 



dn m 
d m {n((j>)) 



n=l 



n=l 



d\P\' 



(6a) 
(6b) 
(6c) 



b1=o+ 



where "f n (4>) represents the step tension of a macro-step, 
which is formed from the merger of n elementary steps, 
and j n {(f>) — 7n(0) + d 2 j n ((j))/d(j) 2 represents the step 
stiffness of the merged step, (n), which describes the 
mean size of the step-droplets, represents the mean num- 
ber of elementary steps in a merged step. Here, (•) rep- 
resents the thermal average. 



B. Monte Carlo calculation 

We used the Monte Carlo method to demonstrate the 
step motion of the nonconserved system [8lj|[46j on a vic- 
inal surface tilted towards the [110] direction (Fig. [2J. 
Then, the chemical potential difference between the am- 
bient phase and the bulk crystal Afi becomes an external 
variable and is the driving force. 

Initially, a macro-step Ar stcp , made of N steps, was 
set almost in the middle of a surface with an area of 
240^2 x 24(h/2 (Fig. 0). The mean surface slope became 
\p\ = A r step /(240v / 2). Periodic boundary conditions were 
imposed in the vertical direction in Fig. [2] For the hori- 
zontal direction in Fig. [5J the left side of the image was 
higher than the right side by N ste p- 

To study the time evolution of the step configuration, 
we adopted a simple Metropolis algorithm. We randomly 
chose a site and allowed its height h(i,j) to increase 
or decrease with equal probability. Then, if the RSOS 
restriction was satisfied, the height was updated by the 
Metropolis algorithm with a probability P described by 



TABLE I: Initial values of external parameters for the Monte 
Carlo calculations. Driving force: A/x/e = ±0.0005. Size: 
240\/2 x 240^2. 



T nknl 

.Label 


Number 
of steps 

A^stcp 


Temperature 
k B T/e 


Surface 
slope 

\p\ a 


Step-step 
attraction 

tint / e 


Al 


1 


0.35 


2.95 x 10" 3 


-0.5 


A24 


24 


0.35 


0.0707 


-0.5 


B24 


24 


0.36 


0.0707 


-0.5 


C24 


24 


0.37 


0.0707 


-0.5 


A240 


240 


0.35 


0.707 


-0.5 


B240 


240 


0.36 


0.707 


-0.5 


C240 


240 


0.37 


0.707 


-0.5 


Al' 6 


1 


0.35 


2.95 x 10 -3 





A24' 6 


24 


0.35 


0.0707 





A240' 6 


240 


0.35 


0.707 






a \p\ = N stcp /(240V2) 
'The original RSOS model 



TABLE II: Monte Carlo results on the surface normal veloc- 
ities and the lag times. Maximum time: 1 x 10 s To- to equals 
1 Monte Carlo step/site. 



Label 


Normal 


Lag time 


Mean size of 




velocity of 




step-droplets 




the surface™ 








v z (xKT 6 di/ro) 


t l /t (xlO 5 ) 


<n> 


Al 


±0.38 ± 0.03 






A24 


±0.63 ±0.15 




8.26 ±0.02 


B24 


±7.30 ±0.08 


8±3 


1.167 ±0.003 


C24 


±7.61 ±0.08 


5±3 


1.102 ±0.002 


A240 


±0.246 ± 0.011 




83.9 ±0.6 


B240 


±22.90 ±0.09 


22 ±4 


5.61 ±0.08 


C240 


±23.87 ±0.16 


7.2 ±0.8 


5.28 ± 0.03 


Al' 


±0.381 ±0.010 






A24' 


±9.07 ±0.06 


0.4 ±0.4 


1.00620 ± 0.00002 


A240' 


±66.2 ±0.3 


0.30 ±0.10 


6.6220 ± 0.0003 



"The velocity with the minus sign corresponds to Afi/e 
-0.0005. 



where AE(i,j) = E(h(i,j) ± 1) - E(h(i,j)) ± A/j, and 
fi = 1/ksT. The energy E(h(i,j j) was calculated using 
the p-RSOS Hamiltonian shown in Eq. Q. The driving 
force of the motion of the crystal surface was designated 
by A/x, where A^i > for growth and A/i < for subli- 
mation. 



p= J 1 (AE(i,j)<0), 
\ exp[-(3AE(i,j)} (AE(i,j)>0), 



The values of the parameters for each simulation are 
shown in Table HI 
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FIG. 2: (Color online) Top views of the vicinal surfaces tilted towards the [110] direction, calculated by the Monte Carlo 
method. The brighter, the higher, with ten gradations. 1 x 10 s MCS/site. A/i/e = 0.0005. e int /e = -0.5 N stcp = 24. Size: 
240^2 x 240^2. (a) A24 (T < 7>, 2 ). (b) B24 (T f>1 < T < 2> >2 ). (c) C24 (T > T f>1 ). (d) A24' (Original RSOS model). 
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1500 
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FIG. 3: (Color online) Time evolution of the mean sur- 
face height under the p-RSOS model. Size: 240^2 x 240^. 
N stcp = 240. A/i/e = ±0.0005. eint/e = -0.5. Full lines: 
A240, k B T/e = 0.35. Broken lines: B240, k B T/e = 0.36. 
Chain lines: C240, k B T/e = 0.37. Two-dot chain lines: 
A240', Original RSOS model (e in t = 0) at k B T/e = 0.35. 



±0.0005, where Af represents the number of lattice 
points. As seen from Fig. [5J the surface motion is linked 
to the step movements. 

Quantitatively, after the time 1 x 10 7 to, where tq 
designates the time for one Monte Carlo step per site 
(MCS/site), the mean surface heights of B24, C24, and 
A24' i ncrea se linearly. We define the surface velocity 
v z as dh(t)/dt. Using the data from the time interval 
2 x 10 7 to to 1 x 10 8 to, we obtained the constant veloci- 
ties by fitting the data to a linear function with the least- 
squares method; these velocities are shown in Table [TT1 
By extrapolating the linear function to the initial height, 
we obtained a finite time T£, which we call the 'lag time'. 
The lag times are also shown in Table HH 

In the case of fceT/e = 0.35 (A24), the center of the 
macro-step hardly moves. 



III. DISCONTINUITY IN THE SURFACE 
TENSION 



C. Normal velocity of the surface 

We present snapshots of the vicinal surface near Tf t i 
and T/,2 at A/i/e = ±0.0005 in Fig. EU From the meth- 
ods of statistical mechanics [3, lilf . we had k B T f l /e = 
0.3610 ± 0.0005 and k B T f ^ 2 /e = 0.3585 ± 0.0007, where 
k B represents the Boltzmann constant. 

In the case of k B T/e = 0.35 (T < T/ )2 ) in Fig. H (a) 
(A24), the merged step hardly moves forward (A/i > 0) 
or backward (A/i < 0). In the case of k B T/e = 0.36 (B24) 
and k B T/e — 0.37 (C24), steps bunch locally, and the two 
cases look similar. However, the mean terrace width of 
B24 looks larger than that of C24. As a comparison, we 
show the case of the original RSOS model (ei n t = 0) in 
Fig. [21(d) (A24'). From the figure, we see that the mean 
terrace width of A24' looks smaller than that of C24. 

In Fig. [3j we show the time evolution of the mean 
surface height h(t) = (l/Af)Yli j M^i) f° r ^/i/e 



A . Calculations from the methods of statistical 
mechanics 

Let us consider a vicinal surface that is so close to 
equilibrium that the diffusivity of atoms does not be- 
come the rate-limiting process[47j]. In that case, the 
interface-limited growth/sublimation becomes the rate- 
limiting process. A detailed understanding of the surface 
thermodynamic quantities are, therefore, essential for un- 
derstanding the behavior of a vicinal surface under small 
\AfJ>\. In this section, we study the discontinuity in the 
vicinal surface free energy f(p) and the surface stiffness 
tensor (f^). 

In order to calculate the surface free energy with the 
methods of statistical mechanics, we add the terms of the 
Andreev fieldQ if = {rj x ,rj y ) as external variables to Eq. 
(0| . The Andreev field tilts the (001) surface to make a 
vicinal surface. The model Hamiltonian given in Eq. Q 
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FIG. 4: For the surface tilted towards the [110] direction calculated by the PWFRG (DMRG) method, (a)-(d) slope dependence 
of the vicinal surface free energy \f{px,p y ) — /(0, 0)]/e, and (e)-(h) surface tension 7surf(p,p)/e . for (a)-(c) and (e)-(g), 
£int/e = —0.5; and . for (d) and (h), ej n t = 0, the original RSOS model. Temperatures: k B T/e = 0.35 for (a), (d), (e), and (h); 
k B T/e = 0.36 for (b) and (f); k B T/e = 0.37 for (c) and (g). Closed squares: (a) and (b), (0,0) and [/(l, 1) - /(0,0)]/e; (e) and 
(f), 7surf(0, 0)/e and 7 su rf(l, l)/ e - Broken lines: (b) and (f), the curves for the metastable states. 




FIG. 5: (Color online) Equilibrium configurations on the p- 
RSOS model. Top views of the vicinal surfaces tilted towards 
the [110] direction calculated by the Monte Carlo method. 
The brighter, the higher, with ten gradations. |p| = 0.707. 
1 x 10 8 MCS/site. Afj, = 0.0 (at equilibrium). e int /e = -0.5. 



step 



= 80. Size: 80^2 x 80^2. (a) k B T/e = 0.35. (b) 



k B T/e = 0.36. (c) k B T/e = 0.37. 



then becomes 

'Hvicinal = 'Hp-RSOS 



-Vy^lHhJ + 1) - h(i,j)] 



(8) 



the Andreev surface free energy as follows: 

f(p) = f(v)+Pn- 



(10) 



Direct calculation of Eq. © is impractical due to 
the complexity of the entropy estimations associated 
with the vast variety of zigzag structures of the sur- 
face steps and with the parallel movements of the steps. 
Fortunately, the density-matrix renormalization group 
(DMRG) method[48[ developed for the ID quantum spin 
system can be used to calculate the partition function 
Eq. © of the p-RSOS model. Using the Suzuki- Trotter 
formula [4{|, a ID quantum spin system can be mapped 
to a transfer matrix[50] of a 2D classical system, such 
as a surface. The transfer-matrix version of the DMRG 
method was developed by Nishino et a/. 51]-[53] fo 



r an 



infinite lattice, and it is called the product wave-function 
renormalization group (PWFRG) method. Since the p- 
RSOS model can be mapped to a transfer matrix, we 
adopted the PWFRG method to calculate Eq. ©. 



B. Equilibrium step configurations 



The partition function Z for the p-RSOS model is given 
by 2 = Y^{h(i .j)} expf-/^ vicinal], where /3 = l/k B T. 
The Andreev surface free energy /(rf)Q is the thermo- 
dynamic potential calculated from the partition function 
Z using 



Pf(v) = ~ J im T7 lnZ ' 
N ->oo J\l 



(9) 



where Af is the number of lattice points on the square 
lattice. The vicinal surface free energy is obtained from 



In Fig. 21 we show the slope dependence of the 
vicinal surface free energy f(p) and the surface ten- 
sion 7 sur f calculated by Eq. (0), Eq. (fl"Uj) . and Eq. 
©, with the PWFRG method. /(0,0) is assumed 
to be e, and 7 SU rf(p,p) represents 7surf(^) with n — 

(-p,-p, i)A/i+i?. 

For k B T/e = 0.35 (Fig. H (a)), only the values of 
the surface tension of the (001) surface and the (111) 
surface were obtained. The surface tension with a mean 
surface slope in the range < \p\ < v2 does not exist 
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for the regular train of steps, because the homogeneous 
surface is thcrmodynamically unstable there. Hence, if 
the mean slope of the vicinal surface has the value of 
< \p\ < the surface becomes a mixture of the (001) 
and (111) surfaces, and its vicinal surface free energy is on 
the tangent line connecting the value of the (001) surface 
with the value of the (111) surface. We demonstrate this 
structure in Fig. [5] (a) using the Monte Carlo method. 

For fc B T/e = 0.36 (Fig. H (b)), the surface tension 
and the vicinal surface free energy increase continuously 
from 7surf(0,0) = /(0, 0) as p increases. Thus, for small 
\p\ or small A^top, a homogeneous structure is expected. 
The homogeneous structure (Fig. H](b)), however, is not 
the same as the regular train of steps seen in the ID FF 
universal system (Fig. [2](d)). The value of (n) is larger 
than one (see Table [TTJ) , which means that step-droplets 
are formed. 

For large 7 su rf(P,p) and f{p,p)/e jump from p t = 
0.349 ± 0.002 to po = 1 at equilibrium. The broken line 
in the figure indicates the metastable state for 0.349 < 
p < 0.501. The tangent line with the end point /(l, 1) 
contacts the f{p,p) curve at f(pt,Pt)- Hence, a vicinal 
surface with mean surface slope p t < p < 1 should be 
formed by the mixture of the surface with the slope p = pt 
and the (111) surface. This structure is demonstrated in 

Fig.HKb). 

For knT/e = 0.37 (Fig. [4] (c)), the surface tension and 
the vicinal surface free energy are continuous for all p. 
Therefore, the mixture of two surfaces does not occur in 
equilibrium (Fig. [5](c)). Step-droplets appear, however, 
due to the sticky character of the steps. 



C. Surface stiffness tensor 

Let us describe a slowly undulating surface by z(x,y), 
where z{x,y) is the surface height at a point (x,y). A 
slowly undulating surface is a surface whose orientation 
varies slowly from place to place around (001). 

Near equilibrium, the time derivative of the height 
of the slowly undulating surface is assumed to equal 
the variational derivative multiplied by a transport co- 
efficient V({p})[13, [EH based on the time-dependent 
Ginsburg-Landau theory. After some manipulation (Ap- 
pendix we have the following equation: 



= v(M)<|^ + E[/^S] 



(11) 



where Q represents the volume of the growth unit in 
the crystal, (/* J ) represents the surface stiffness tensor 

d 2 f{p)/dp i dp j \ i ; = p c (p= (p x ,p y ) = (pi,P2))[H|, and z^f 
represents d 2 z jdx 1 dx^ (x = {x,y) = (x 1 ^ 2 )). Here, p e is 
the surface slope at equilibrium, and df(p)/dpiy = .p — 
(< = {1,2}). 

The second term on the right-hand side of Eq. pT|) 
expresses the x-y projected Gibbs- Thomson effect. From 



the condition v z = 0, we obtain another expression of the 
ECS, except for the facets. 

For T > T L2 , ksT/e = 0.36, and k B T/e = 0.37, the 
surface free energy f(p) increases monotonically as \p\ 
increases around \p\ — (Fig. [4] (b) and (c)). We can 
then calculate the surface stiffness tensor (/ iJ ) in Eq. 
(fTTj) explicitly in the limit of \p\ — > 0. 

Keeping the non-GMPT \p\ expanded form of the vic- 
inal surface free energy (Eq. (|3])) in mind, we describe 
the vicinal surface free energy as follows: 

f(f) = /(O) + 7W>)f + + o(Ip1 c (12) 

where d represents the unit height of a single step (d — 
1), j(4>) represents the step tension of a single step, and 
B^{4>) represents the coefficient of |p|^/c?^. 

Adopting Eq. (TT2")) . and after some calculations, we 
obtain the expressions of the surface stiffness tensor 
(/ IJ ')[56| in the limit of |$ as follows: 



f U = 



|p1 C " 2 (ii 



t2 sin 



r = 



m . 2 

. -: sin 
\P\ 

-2t 3 sin0cos0) +0{\p\ c - 1 ), 

7W_2^ , |flC-2(f 4 + t2aill 2^ 

6) + o(bf- 1 ), 

|pf- 2 Ml-2sin 2 



+2*3 sin <fi cos 



[p\ 



■ sin cos ( 



+t 2 sin^cos^] +0(\p\ c - 1 ) 



(13a) 



(13b) 



(13c) 



where 7(0) (j((f>) — "f{cj)) + d 2 j((j))/d(j) 2 ) represents the 
step stiffness of a single step, and t±, t 2l £3, and t± are as 
follows: 



(14) 



h 


= c(C-i)B f (0), 


t 2 


= C(C- 2)^(0) 


h 


= (C-i)B' c (4>), 


u 






= d 2 B ( {(j))/d(j) 2 . 



The principal values of the become 
fn = bf-^+Odpf" 1 ), 

ft = W- + m c - 2 u+om c - 1 )- 
\p\ 

Then, we obtain dct(f^) as follows, 

det(,P) = h\p\<- 3 m + 0(\p\^ 2 ) 

= \p\^ 3 ac^i)B C m(<p) 

(\P\ -+ 0). 



(15a) 
(15b) 



(16) 



Physically, f n represents the surface stiffness against a 
bending stress that is normal to the facet edge, and ft 
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represents the surface stiffness against a bending stress 
that is tangent to the facet edge. 

Using /„ and ft , we have an expression for the normal 
surface velocity, as follows: 



where zffl and z^' are the eigenvalues of {d 2 z / dx % dx 3 ) 
(i,j — 1,2). According to recent dev elop ments in the 
study of nonequilibrium bunched stepsj57|-[59j . the pro- 
file of a bunched step is related to the force range of the 
effective step-step interactions on the nonequilibrium vic- 
inal surface. 



V(p) 



,(2) 



(17) 



At equilibrium, where v z 
obtained explicitly as 



0, the 4 2) 



and Zi are 



z n % = -[A/V(2fi)]b1 2 -V[C(C- i)B C m, 



(2) 
z t,eq 



= -[A/V(2fi)]|p1/7(0) 



(18) 



in the limit of \p\ — » 0. Therefore, the Gaussian curvature 
K near the (001) facet edge on the ECS is obtained as 
follows: 



A 2 



A 2 b1 



3-C 



K ,9 2 det(/y) ~ C(C-1)B{(^)' (b1 ^ 0) 

(19) 

where 5 = 1 + \p\ 2 , and A = A/u/(2f2). 

Therefore, for Tf^ < T < T/ i, applying the form of 
vicinal surface free energy (Eq. (J3J) to Eq. (TT2"j) and 
Eq. (|15ap ~Eq. ([T9"l) . we have the principal values and 
the determinant of the stiffness tensor and the Gaussian 
curvature as follows: 

f n = 2A eB (4>), f t =^(<f>)/\p\ (20a) 
det(r')= 2A eS (4>)^ 1 (4>)/\p\ (20b) 
K = \ 2 \p\/{2A cS (m<P)} (20c) 

where C = 2[Hj]. 

For T fA < T « T R , where C = 30 and T R is 
the roughening transition temperature of (001) surface 
( WA|) . we have 

/„ = 6B eff (^M / t = 7 i(0)/b1 (21a) 
det(/«)=6S eff (^) 7l (0)/N (21b) 
K = A 2 /[6B cff (0) 7 (0)]. (21c) 

Due to the step-droplets, B e g{(j)) (Eq. (l5bf ) has a differ- 
ent value from the value of the ID FF[15] in the following: 



Bi(^) = (k B Tn) 2 /[6^)} 



(22) 



For T/,i << T < Tr, the p-RSOS system converges 
to the ID FF system; i.e., i? c ff(</>) converges to B 1 ((f)). 
Hence, det(.P 3 ') and if converge to the GMPT universal 
value 0] of 

det(f ij ) = {k B Tnf and K = \ 2 /{k B Tir) 2 , (23) 



TABLE III: Surface stiffness and Gaussian curvature. 



Lable 




Pft b 




iv/(/3A) 2d 


B24 


1.24 x 1(T 2 


19.4 


0.241 


4.10 


C24 


0.119 


19.3 


2.30 


0.431 


A24' 


0.507 


19.5 




l/Gpr) 2 ' 








« 9.870 


w 0.1003 



»j8 = l/fc B T. Eq. ifTSat 

6 Eq. jug 

c Eq. JTS} 
d Eq. JT5} 

e GMPT (ID FF) universal value, ref.[H 



1 + |p| 2 . |p| = 0.0707. 



respectively. 

In Table IIII1 we show the values of the surface thermo- 
dynamic quantities for k B T/e — 0.36 (B24) and k B T /e = 
0.37 (C24). As a comparison, we show the values of the 
original RSOS model at k B T/e = 0.35 (A24'). In the 
process of the estimation of f n and f t , we used values 
calculated by the PWFRG method for f3A ctf (ir/A) and 
/3Beff(7r/4)0; i.e., pA cS (n/4) = (6.22 ± 0.06) x lO" 3 
for B24 and /3B eS (ir/4:) = 0.281 ± 0.008 for C24. We 
also approximated /3 7 i(7r/4) by the values of the inter- 
face stiffness of the 2D nn Ising model[60-62]; i.e., 1.377, 
1.371, and 1.364 for A24', B24, and C24, respectively. 

As seen from Table IIII1 even though f t is kept al- 
most constant, /„ changes drastically as the tempera- 
ture changes. This means that the step-droplets soften 
the surface against the bending force normal to the mean 
running direction of the steps. In fact, the Gaussian cur- 
vature at the (001) facet edge increases. 

Note that p on the ECS (z — z(x, y)) depends on (x, y). 
Taking the s-axis normal to the edge of the (001) facet 
at a point (x c , 0), 



P\[z(x,0) - z(x c ,0)} = -An[/3\(x - x c )} 6 



(24) 



where j3 = l/k B T, 9 n is the normal shape exponent, and 
A n is the normal amplitude. Using Eq. (fl~2"|) . 9 n and A n 
arc then described as follows! 191: 



c 



c-i 

X(x - x c ) 



l 



k B T 



(25a) 
(25b) 



IV. 



INHIBITION OF STEP MOTION DUE TO 
'STEP FACETING' 



A. Roughness on the side surface of a merged step 

In the section, we explain the inhibition of the macro- 
step motion in the vicinal surface at sufficiently low tem- 
peratures T < Tf } 2- 
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For T < T/,2, the inhomogeneous vicinal surface, a 
mixture of the (001) and (111) surfaces (Fig. [5] (a)), is 
realized due to discontinuities in the surface tension [19j. 
Since both the (001) surface and the (111) surface are 
smooth, a macro-step in the vicinal surface becomes a 
faceted step ('step faceting' [43] ). Hence, the squared sur- 
face widths for the (001) and (111) surfaces should both 
be finite. 




0.8 0.9 1 

(xlO 8 ) 

MCS/site 

FIG. 6: Intermittent sublimation of the surface. A24. 
Afi/e = -0.0005. e int /e = -0.5. k B T/e = 0.35. N stap = 240. 
Size: 240\/2 x 240\/2. 

The squared surface width of a vicinal surface W(n) 2 
is defined as follows [56j: 

W{nf = {[z(x) - (z{x))?)/9, (26) 

where (•) represents the thermal average and n — 
{—pi, — P2j l)/<7 represents the normal unit vector of the 
tilted surface. From the result of Ref.[56|, W(n) 2 is 
connected to the surface stiffness tensor in the limit of 
L — > oo, where L is the linear size of the surface, as 
follows: 

W(nf = kBT • In L. (27) 

2ng^/det(fv) 

Since W(n) 2 should be finite for a smooth surface, 
det(/ y ) should be divergent in the order of (InL) 2 in the 
thermodynamic limit. Hence, slight deformations from 
the flat surfaces of (001) or (111) are pulled back to the 
flat surfaces because of the strong Gibbs-Thomson effect. 

(2) (2) 

Further, if we consider z n = and z\ ' =0 on 
the macroscopic vicinal surface, the transport coefficient 
V(p) reduces to zero because the kink density on the side 
of the macro-step converges to zero in the thermody- 
namic limit. Therefore, the continuous motion of the 
surface described by Eq. (fTTj) or Eq. p7|) is inhibited in 
this temperature region. 

Instead, as shown in Fig. [51 the intermittent motion of 
the surface occurs in a longer time scale. The mechanisms 
that cause the intermittent surface motion are considered 
to be 2D nucleation j4] around the intersection line of two 
surfaces, a terrace surface and the side surface of a faceted 



step. Near equilibrium, the nucleation rates at about the 
center of the terrace and at about the center of the side 
surface are small, and the nucleation rate around the 
intersection line of the surfaces becomes larger because 
of the special geometrical arrangement. 



B. Self-pinning of steps 




FIG. 7: (Color online) Self-pinning of steps on the p- 
RSOS model. Top views of the vicinal surfaces tilted to- 
wards the [110] direction calculated by the Monte Carlo 
method. The brighter, the higher, with ten gradations. 
Afi/e = 0.1. e int /e = -0.5. k B T/e = 0.1. iV s t cp = 24. Size: 
240v / 2x240v / 2- (a) Initial configuration, (b) 1000 MCS/site. 
(c) 2000 MCS/site. (d) 3000 MCS/site. (e) 5000 MCS/site. 
(f) 1 x 10 4 MCS/site. 
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(a) 




1 


in 

91 


I 
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1 





FIG. 8: (Color online) Step growth on the p-RSOS model. 
Top views of the vicinal surfaces tilted towards the [110] direc- 
tion demonstrated by the Monte Carlo method. The brighter, 
the higher, with ten gradations. A/i/e = 0.35. eint/e = —0.5 
fc B T/e = 0.1. TVstcp = 24. Size: 240^ x 240^- (a) Initial 
configuration, (b) 1000 MCS/site. (c) 3000 MCS/site. (d) 
5000 MCS/site. (e) 9000 MCS/site. (f) 1.3 x 10 4 MCS/site. 
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At sufficiently low temperature, the step movements 
of the macro-steps are almost inhibited for small A/i be- 
cause of step faceting. In Fig. [Jj we show the time 
evolution of the vicinal surface at k B T /e = 0.1 for the 
p-RSOS model with ej nt /e = —0.5, starting from the 
regular train of elementary steps (Fig. [7] (a)). In our 
previous paper [l8j. we showed that for A/z = 0, step zip- 
ping, the phenomena of successive sticking of adjacent 
steps, much like a slide fastener in clothes, occurs at the 
collision points of adjacent steps. 

In the case of A/i/e = 0.1, step- flow growth occurs in 
the early stages (Fig. [7] (b)) of surface growth. When 
an elementary step collides with an adjacent elementary 
step, the steps merge to form a double step due to the 
sticky character of steps (Fig. [7] (b)). Since the step- 
droplet now has a velocity lower than that of the ele- 
mentary steps, the steps from behind catch up with the 
new step-droplet and merge to form a larger step-droplet 
(Fig. 0(c)). Since the larger step-droplet now has even 
lower velocity, step movements become pinned by step- 
droplets, as in a traffic jam (Fig. [7J(c)-(e)). The mobility 
of the merged steps is so slow that the vicinal surface is 
almost quenched when elementary steps disappear (Fig. 
EJe),(f)). This inhibition of the step motion occurs when 
A/x/e < 0.3 and around A/i/e = 0.3 if k B T/e = 0.1. 

At A/i/e = 0.35, elementary steps begin to separate 
successively from the lower side of the macro-steps (Fig. 
[8]). This separation of elementary steps looks like step 
nucleation. As mentioned at the end of the previous sub- 
section, step nucleation starts from the 2D nucleation 
around the intersection line between a terrace surface and 
a side surface of a macro-step. In addition, the growing 
steps in Fig. [5] exhibit distinctive wavy shapes due to 
the sticky character of the steps. The frequency of the 
separation of steps is not simply defined, and it may be 
considered to be related to 'kinetic roughening ' [3, 1(33 . [64 1 
on the surface. 

This separation of steps occurs more frequently when 
A/i/e > 0.35, so we designate A/i* /e = 0.35 the threshold 
driving force for the separation of a step at k B T/e = 0.1. 
When the temperature increases within the range of 
T < 7/ 2, the value of the threshold force for the sep- 
aration of step A/i* becomes smaller. This is because a 
decrease of the step tension for an elementary step, which 
is caused by the entropy of the zigzag structure in an el- 
ementary step, induces more frequent 2D nucleation at 
the intersection line. 

In the case of sublimation (A/i < 0), an elementary 
step appears from the upper side of the macro-step, and 
the step goes backward. This separation of elementary 
steps also looks like 'step nucleation', and the steps that 
move backward exhibit distinctive wavy shapes, which 
are almost mirror-symmetric to the shapes observed for 
the growing steps. 



TABLE IV: Surface roughness and surface velocities. A/i/e 
±0.0005. 



Label W{n) 2 /\nL a 


On' 
Hk,x 


c 

Pk,2 


d 

Pk 


v e (x 10~ 6 rfi /m) 


Al 





0.116 




0.116 


0.346 


B24 


0.322 


0.118 


0.0598 


0.108 


7.49 


C24 


0.104 


0.119 


0.0622 


0.113 


7.63 


Al' 




0.131 




0.131 


0.388 




w 0.05066 










A24' 


1/(2tt 2 3) 


0.131 


0.226 


0.131 


9.36 




w 0.05041 











"Eq. j27} 

b Eq. 03} 
c Eq. jUJ 
d Eq. Eg} 

e Eq. 02} and Eq. J38} 

fg = l + \p\ 2 



V. SLOWING DOWN OF THE STEP MOTION 
DUE TO STEP-DROPLETS 

A. Roughness of the vicinal surface 

In this section, we study the step motion in the case 
of 7/ 2 < T ~ Tf t i. In this temperature range, the sur- 
face tension and the vicinal surface free energy f(p) are 
discontinuous at the (111) surface but continuous around 
\p\ = 0. Hence, Eq. (fTTj) is applicable, and it describes 
the motion of the vicinal surface around \p\ = 0. Also, 
for the mean flat surface, the velocity of the surface v z is 
expressed using Eq. (ITTT) as follows: 



(28) 



According to simple linear response theory [65l|. the trans- 
port coefficient V(p) should be proportional to the 
squared fluctuation width W 2 (n) of the vicinal surface 
due to the fluctuation-dissipation theorem (Appendix [Bl 
Eq. (|B6[) h In this subsection, therefore, we study the 
roughness of the vicinal surface. 

Let us first evaluate the roughening transition tem- 
perature T R of the (001) surface. By using the PWFRG 
method, we calculated T R from the universal relationship 
K R /\ 2 = 4/(k B T R 7r) 2 , where K R is the GMPT Gaus- 
sian curvature on the ECS at T = T R 0, [H] . The ob- 
tained T R are k B T R /e = 1.505 ± 0.008 for e int /e = -0.5 



and k B T 



RSOS 



1.584 ± 0.006 for e int = 0. The 



value of Tn SOS is consistent with the value obtained by 
den Nijs(66|. The result that T R < T^ sos means that 
the step-step attraction slightly roughens the (001) sur- 
face. The local bonds e- mt are considered to stabilize the 
'blobs' [671] in the zigzag structure of a single step. These 
blobs enhance the deformations of a step by thermal fluc- 
tuation, which decreases the step free energy per length. 

The roughness of the vicinal surface is measured by the 
squared surface width W(n) 2 . Hence, by substituting Eq. 
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(fT6|) into Eq. (|27|). we have 



W(n) 5 



27n,VC(C - 1) WHWO 



•InL. 



(29) 



Using the values of det(/ IJ ) in Table lllll we calculated the 
corresponding values of W{n) 2 / InL and present them in 
Table ED 

For T /i2 < T < T f>1 , we can see from Table [IV] 
that the step-droplets roughen the vicinal surface. For 
< T << Tn, the vicinal surface is still roughened 
by the remaining step-droplets. As the temperature in- 
creases from 2/i, the roughness of the vicinal surface 
W(n) 2 /\nL (\p\ -» 0) rapidly converges to the GMPT 
value 1/(2tt 2 ). 

In summary, for the cases demonstrated in £|II C[ we 
have the following relationship from the values in Table 

El 



W 2 {n)\ f= ^ & >W 2 {n)\ f= 



0.37 



> W 2 (n) 



2/-^RSOS| 



lT=0.37 



W 2 (n) 



RSOSi 



lT=0.35' 



(30) 



where T represents the reduced temperature k^T /e. 




□ □SB 



4exp[-|3e] 4exp[-|3s] 4exp[-|3s] 2exp[-2[Ss] 4exp[-|3(2E+E. )] 




statistical 
weight 

number 
of kinks 



FIG. 9: (Color online) Conceptual diagrams of kink config- 
urations for an elementary step. A/i > 0. Top figure: an 
example of an elementary step (top view). The lattice sites 
are divided into two sub-lattices, distinguished by stripes and 
solid. Darkly shaded areas describe the lower terrace. 'K' 
represents a kink site. 



B. Step smoothing 



From Eq. (f3T)]) . referring to the expression of V(p) 

> 



(Eq. (|B6|) h we must have the relationship v z \f =0 



v z\±=0.37 V 



RSOSl 



3(i 



IT=0.35 



for the cases demonstrated in 
j jll CI However, the calculated order of surface mobility 
is clearly seen from Fig. [3] and from Table [TTI as 



^lT=0.36 



< v z\t 



T=0.37 



< v RSOS \ ' 

^ u z lT=0.35' 



(31) 



Therefore, even near equilibrium, the simple linear re- 
sponse theory is not applicable to describe the transport 
coefficient for the vicinal surface with sticky steps around 
T~T fjl . 

It is the formation of the step-droplets that prevents 
the application of the simple linear response theory to the 
transport coefficient. Though the step-droplets roughen 
the vicinal surface, they also diminish the number of 
kinks on the side of the merged steps. At temperatures 
sufficiently lower than Tr, the step motion is governed 
by the total in/out flow of materials at the kink sites. 
We will use 'step smoothing' to refer to this reduction in 
the number of kinks. 

We note that the transport coefficient is described as 
follows: 



V(p) = Pk N [lw] N step n 2 /(Sk B T), 



(32) 



where pk represents a kink density, iV[j 10 ] represents the 
linear number of the unit cell in the direction of [110], and 
S = 240 x 240 x 2 represents the size of the simulated area. 
We show a 'kink' in the case of the growth (A/i > 0) of 
a single elementary step in Fig. |H] The kink density for 



a single elementary step is approximated based on the 
19- vertex model[43] as follows (Fig. 



Pk,i 



(1 + 2e-^) 2 
[6 + c-/ 3e + 2e-' 3 ( e + e "«)]2- 



(33) 



Using the expressions of pk,i and V(p), we obtained v z 
for an elementary step, as shown in Table HVl for the data 
of Al and Al'. By comparing these values with the ones 
in Table [TTI it is seen that the values obtained by Eq. (|3"2"]l 
and Eq. (I33[) are in agreement with the values obtained 
by the Monte Carlo method without fitting parameters. 

When \p\ increases, that is, 7V s t op increases, the step- 
droplets are created by a local merging among steps. 
Since (n), which represents the mean size of the step- 
droplets (boson n-mers), is close to 1[19| for the surface 
with iV s tep = 24 (Table |TT]), the double step is statisti- 
cally dominant among multiply merged steps. Hence, we 
consider the mean kink density as follows: 



p k = {2- (n))pk,i + ((n) - l)pfc,2, 



(34) 



where represents the kink density of the double step. 
In equation (fM)) . we assumed (n) s» 1-(1— p)+2-p, where 
p represents the fraction of steps that are double steps. 

Since the configuration of a double step (Fig. [TUJ) is 
approximately described by the triple vertices, pk,2 is es- 
timated as follows: 

[2 + 2e~ /3 ( e+£int ) -I- e -£(2e+2£i„t)] e /3eint 

PK2 ~ 3[3 + 2e-^+^ + e-^ 2 ^ 2 ^)] ' ^ 

In Table IIV1 we show the values of Pk.i, Pk,2, & n d Pk 
calculated by using Eq. ([55]). Eq. ([S5])'. and Eq. (jMl) . 
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statistical 
weight exp[-(3(4e+E int )] 

number 
of links 



LT 

exp[-P(4s+e int )] 



1 



1 



exp[-P(4 £ +£ int )] 





ft 




exp[-P(5E+2E int )] exp[-P(5E+2s int )] exp[-P(6£+3e int )] 
1 1 1 



FIG. 10: Conceptual diagrams of kink configurations for a 
double step. A/i > 0. Top figure: an example of a double step 
(top view) . More darkly shaded areas describe lower terraces: 
white > gray > gray with stripes > gray with hatching > 
black. 'K' represents a kink site. 



respectively. Here, we used the value of (n) in Table [TT1 
Comparing the values of v z in Table ITVl obtained by Eqs. 
(H?2"1) - (|35l) with the ones obtained by the Monte Carlo cal- 
culations ( mi Cp . we see that they are in close agreement 
without fitting parameters. 



VI. SUMMARY AND DISCUSSION 

We studied the step motions in a vicinal surface near 
equilibrium using a lattice model with sticky steps. We 
adopted a restricted solid-on-solid model with a point- 
contact- type step-step attraction (Eq. (JU), e- mi < 0, the 
p-RSOS model). We showed that the point-contact-type 
step-step attraction caused discontinuities in the surface 
tension 7 sur f(n) (Eq. @) and in the vicinal surface free 
energy f(p) at low temperatures (Fig. H](a)). Due to the 
discontinuities, 'step faceting' [42] occurred on the macro- 
steps in the vicinal surface tilted from the (001) surface 
towards the [110] direction for T < T f>2 (Fig. [5] (a)). 
The continuous motion of the macro-steps was inhibited 
under a small driving force A/i by the step faceting (Fig. 
[5] (a), Fig. [3]). Instead, the intermittent motion of the 
vicinal surface (Fig. [SI Fig. [5J occurred by way of 2D nu- 
cleation around the intersection line between the terrace 



surface and the side surface of the macro-step. Starting 
from the regular train of steps as the initial configuration, 
pinning of steps was demonstrated to take place from the 
collision point of the adjacent elementary steps (Fig. [7]) 
without impurities or defects. We term this phenomenon 
the 'self-pinning' of steps. 

For < T < Tf i, the surface tension and the vic- 
inal surface free energy were continuous for the surface 
slope \p\ ~ 0, but they were still discontinuous around 
the (111) surface (Fig. H](b)). This discontinuity around 
the (111) surface led to the formation of 'step-droplets' 
(boson n-mers), locally merged steps. The step-droplets 
roughen the vicinal surface, and we showed this by cal- 
culating the squared surface width W 2 (n) fflVAl Table 
irV|) . using the relationship between the squared surface 
width and the determinant of the surface stiffness ten- 
sor det(/ 4i ) ( jnrgl Table HB. The step velocity, which 
is expected to be larger than that of the original RSOS 
model, was lower than the step velocity of the original 
RSOS model at the same temperature (Fig. [3l Table ITT1) . 
This is because the step-droplets diminish the kink den- 
sity of the vicinal surface (the 'step smoothing', WBI) . 
We estimated the transport coefficient (Eq. ([3"2"j) and 
([34|)) using the kink density of an elementary step pk,i 
(Eq. (J22])) and the kink density of a double step pk,2 
(Eq. (1551) ). Using the equations of the kink densities, we 
reproduced the surface velocities obtained by the Monte 
Carlo method ( mi C[ Table IIV[) without fitting parame- 
ters. 

The time lag of the surface motion for an abrupt re- 
verse of the driving force around equilibrium is one of the 
methods of detection for the step-droplets that are cre- 
ated by the discontinuity in the surface tension. When 
A/i > (growth) at low temperature, elementary steps 
separate successively from the lower side of the step edge 
on a macro-step (Fig. |8](b), (c)). The elementary steps 
in the terrace of the upper side of the macro-step catch 
up and merge with the macro-step (Fig. [5] (d), (e)). 
Hence, the profile of a macro-step is asymmetric with 
respect to the upper side and the lower side of the step 
edge. Then, when the driving force is reversed abruptly 
to A/i < (sublimation), the movement of the elemen- 
tary steps reverses; for the macro-step, however, the pro- 
file of the step edges changes first. In this way, a time 
lag of the surface motion occurs. If step droplets do not 
exist, as in the original RSOS model, the bunched steps 
dissolve rapidly because the edges of the side surface of 
the bunched step are rough near equilibrium. There is 
no nucleation barrier for the separation of steps from the 
bunched steps. Therefore, the inhibition of the macro- 
step motion against the alternative change of A/i around 
A/i = becomes the evidence of the discontinuity in the 
surface tension. 

It should be noted that under equilibrium, step merg- 
ing or step faceting do not always occur for sticky steps. 
In order for step merging to occur under equilibrium, 
discontinuities in the surface tension and in the vici- 
nal surface free energy are essential. For temperatures 



12 



higher than in the case of the p-RSOS model, dis- 
continuities disappear because the entropic repulsion be- 
tween the steps overwhelms the short-range step-step 
attraction. Step faceting then disappears. The step- 
droplets, on the other hand, remain for temperatures 
slightly above T/,1, due to the finite character of the step- 
droplets. Therefore, a slowing down of the step velocity 
occurs when |A/z| is small. 

Recently, an anomaly together with strong anisotropy 
in the stiffness of vicinal surfaces was observed [68( for the 
surface around the (0001) facet of a 4 He crystal particle 
at low temperature. This anomaly in the stiffness of the 
4 He crystal particle may be explained by anomalous be- 
havior of the surface stiffness tensor that originated from 
a discontinuity in the surface tension ( mil Q i flVI) . Based 
on this present work and our previous work, we consider 
that short-range step-step attraction with a microscopic 
origin might exist on the surface of a 4 He crystal. A 
hydrodynamic interaction of steps [6^] may also be con- 
ceivable for He. Since the hydrodynamic interaction of 
steps is long range, it causes a discontinuity in the sur- 
face tension of the vicinal surface tilted in all directions 
around the (0001) facet at a temperature lower than the 
specific temperature Tf. In the case of the point-contact- 
type step-step attraction, however, the attraction causes 
a discontinuity in the surface tension only around the 
vicinal surface tilted towards a special direction where 
the step-droplets are formed. In the p-RSOS model, the 
vicinal surface tilted towards the [100] direction from the 
(001) surface shows typical GMPT behavior[19j. We ex- 
pect further relevant experimental studies in the future. 

In the present Monte Carlo simulation, only non- 
conserved attachments and detachments of atoms [46j 
were taken into consideration in order to show clearly 
the effect of the point-contact-type step-step attraction. 
That is, other effects that occur on a real surface were ig- 
nored, such as surface diffusion 22 1, electromigration[27|- 
[29l |. the Schwoebel effect pit [30| . the Shockwave effect |3lj|. 
impurity effects l32|-l34l. s train effects [35[-[4l|, and the ef- 
fect of surfactants |70l |71|. On real surfaces, these effects 
exist together with discontinuities in the surface tension. 
For example, for a system with short-range attraction 
and long-range repulsive step-step interaction, such as 
an elastic interaction, Tersoff et al. [361 ] , Liu et al. 37] , and 
Shenoy et aZ.[38[ showed that a regular array of merged 
steps of size n appears on the vicinal surface, where n 
depends on the strength of the step-step attraction and 
on the strength of the elastic repulsion. 

The vicinal surface of Cu (lln) (n= 5,7,9)[z2 may 
be an example of a surface with a discontinuity in the 
surface tension together with an Ehrlich-Schwobel (ES) 
barrier. The STM image (Fig. 2 (a)) of the Cu (lln) 
vicinal surface in the paper of Neel et al. looks similar 
to the pattern shown in Fig. [7] of the present work. The 
dark spots in their STM image, which look like holes or 
trenches, seem to pin the macro-steps. An explanation 
for the pinned phenomena was not given in their paper. 
We consider that the trenches may be formed by the self- 



pinning caused by the step faceting in the early stage of 
step bunching. The authors said that the step bunch- 
ing had a kinematical cause because the regular train of 
steps was observed at about 700 K. If the characteris- 
tic temperature Tf t i is lower than 700 K, however, the 
macro-steps caused by the sticky character of the steps 
dissolve at T ~ 700 K. 

In our previous papers [z!, 174} . we presented lattice 
models to describe the vicinal surface with adsorption. 
We showed that the discontinuity in the surface tension 
and the vicinal surface free energy is induced by adsor- 
bates. Hence, the inhibition of step motion around equi- 
librium, as presented in this work, is expected to occur 
on a vicinal surface with adsorbates. In fact, the shape 
shown in Fig. [5] and intermittent growth, as is one shown 
in Fig. [BJ are similar to those seen in the observations of a 
Au/Si(001) surfacefnl and an 0/Ni(977) surface^. In 
the case of the Au/Si(001) surface, gold-induced faceting 
occurs, and the origin of the faceting is thought to be 
discontinuities in the surface tension [76j. 

In the case of the 0/Ni(977) surface, oxygen- induced 
step merging has been observed, and it has also been ob- 
served on the surfaces of several other metals [77}. So far, 
step bunching, or step merging, has been considered to 
be a dynamical phenomenon that takesplace when the 
surface is far from equilibrium [5!|, [27]- 34], [78]. In their 
studies of step bunching resulting from various causes 
when the system is far from equilibrium, the surface free 
energy and the step free energy are assumed to have the 
GMPT universal form (Eq. flU). We think that the 
oxygen-induced step merging that occurs near equilib- 
rium may be explained by the discontinuity in the sur- 
face tension that is induced by adsorbates. For a full 
understanding of step dynamics on the vicinal surface in 
combination with step-step attraction and other surface 
effects, further study is required. 



VII. CONCLUSION 

Point-contact-type step-step attraction causes discon- 
tinuities in the surface tension j SUI f(n) and in the vicinal 
surface free energy f(p) around the (001) and (111) sur- 
faces at low temperature. Due to these discontinuities, 
'step faceting' occurs, which is where the side surface of 
a macro-step becomes smooth. Step faceting inhibits the 
continuous motion of macro-steps under a small driving 
force Afi. Step faceting also induces intermittent motion 
of the surface in the long-time behavior, and it induces 
the pinning of steps even if there are no impurities, ad- 
sorbates, or defects on the surface. 

For temperatures such that the discontinuity in the 
surface tension occurs only around the (111) surface, the 
slowing down of step movements occurs due to 'step- 
droplets' (boson n-mers with finite lifetimes), which are 
locally merged steps. The step-droplets roughen the vic- 
inal surface, but they diminish the kink density of the 
vicinal surface (step smoothing). 
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From Eq. (|A1|) and Eq. (|A7[) . therefore, we have Eq. 

CD- 

Appendix B: Transport coefficient based on the 
simple linear response theory 



Let us consider the Langevin equation for A/i = as 
follows 79): 



Appendix A: Equation of surface motion near 
equilibrium 

In this appendix, we derive an equation for the surface 
motion near equilibrium. 

Based on the linear response theory near 
equilibrium[65], the time derivative of the height 
change of the vicinal surface is assumed to equal 
the variational derivative multiplied by a transport 
coefficient. Let z'(x,y) = z(x.y) + 5£5(x,y), where 
5{x — Xq, y — uq) is the delta function. Then, we write a 
kinetic equation for the height change of the surface as 
follows 541: 



St 



-v(W) 



SG 



(Al) 



where V({p\) is a transport coefficient that depends on 
the surface gradient and G is the total free energy of the 
system. The functional G for z'(x,y) is described as 



G = 



dxdy[m-^z\x,y)}, (A2) 



where p = (pijPa) = (dz'/dx,dz'/dy), and f{p) repre- 
sents the surface free energy per x-y projected area (vic- 
inal surface free energy). f(p) depends on the surface 
tension 7 surf (n) as f(p) = 7 sux f (ft) ^Jl + p\ + p|j where ft 
represents the unit normal vector at the surface. 

Considering SG = 0, where the free energy becomes 
minimum, at equilibrium, we have 



j^s/(p) + a/i =o om 

dx 1 - dpi Q dp 



(A3) 



where x = (a; 1 , a; 2 ) = (x,y). From the solution of Eq. 
[3|, we obtain 



ff= -[An/(2fl)]S, f(ff) = /(j% =#e -ff-pe. (A4) 

Then the relation Eq. (|A3|) is inverted as 

Pe = -dfffi/drf. (A5) 

Expanding G around the equilibrium p = p e , we write 
SG as 



5G = 



dxdy 



d 2 f(p) 



dpidpj 



Axdy 



d 2 f(p) 



dpidpj 



P=Pe 

d 2 z 



SpiSp-j - 



(A6) 
(A7) 



^ = uV 2 h + R(t), 



(Bl) 



where we assume R(t) is a Gaussian white noise such as 
(R(t)R(t')) = 2D h S(t - t'). Then, the distribution func- 
tion P(h, t) of the Fokker-Plank equation is proportional 
to 



exp 



d 2 x- 



2D, 



-\{p-fe 



(B2) 



Comparing Eq. (|B2[) with Eq. (|A6|) at Ap, — 0, where 
Eq. (|A6I) should be the Hamiltonian of the Fokker-Plank 
equation, we see that /W) /k^T in Eq. (|A6| corresponds 
to v/2Dh in Eq. (|B2|) . Hence, we have 



V(p) = 2D h /k B T. 



Since Rj cc (pj — Pj >e ), we have 



(RiRj) = const, x k B T(f (af3) -%/L 2 



from which follows 



(B3) 



(B4) 



2 ^ ^ 2 g 2 (W{ftff 



det(R l Rj) = const, x 



L 4 (lnL) : 



(B5) 



Therefore, we approximate V(p) as follows: 



V{p) w const, x 



47r 3 VF 2 (n) 
/s B rL 2 InL' 



(B6) 



Appendix C: Kink sites for sublimation 

In the case of sublimation (A/x < 0) at low temper- 
ature, 'atoms' (elementary cubes in the RSOS model) 
detach mainly from the kink sites shown in Fig. 1111 
Equations for the kink densities for an elementary step 
and for a double step are the same as Eq. (|3"3"]l and (|35l) . 
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FIG. 11: Conceptual diagrams of the kink configuration for 
A/i < (sublimation). Top line: kink configurations of an 
elementary step. Second and third lines: kink configurations 
of a double step. More darkly shaded areas describe the lower 
terraces: white > gray > gray with stripes > gray with hatch- 
ing > black. 'K' represents a kink site. 
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